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Abstract 

We explore how inherent flexibility of a protein molecule influences the mechanism controlling 
the kinetics of allosteric transitions using a variational model inspired from work in protein folding. 
The striking differences in the predicted transition mechanism for the opening of the two domains 
of calmodulin (CaM) emphasizes that inherent flexibility is key to understanding the complex 



r^ [ conformational changes that occur in proteins. In particular, the C-terminal domain of CaM 

^^' (cCaM) which is inherently less flexible than its N-terminal domain (nCaM) reveals "cracking" 

or local partial unfolding during the open/closed transition. This result is in harmony with the 

>, 

0^ ! picture that cracking relieves local stresses due to conformational deformations of a sufficiently 
ON 

^«0 ' rigid protein. We also compare the conformational transition in a recently studied "even-odd" 

t~^ ■ paired fragment of CaM. Our results rationalize the different relative binding affinities of the EF- 

o. 

OO , hands in the engineered fragment compared to the intact "odd-even" paired EF-hands (nCaM and 



cCaM) in terms of changes in flexibility along the transition route. Aside from elucidating general 



^ ■ theoretical ideas about the cracking mechanism, these studies also emphasize how the remarkable 

intrinsic plasticity of CaM underlies conformational dynamics essential for its diverse functions. 



INTRODUCTION 

To understand protein function, it is often essential to characterize large conformational 
changes that occur upon ligand binding. Within the population-shift mechanism of allosteric 
conformational change,- the bound complex is formed when a ligand selects and stabilizes 
a weakly populated conformational ensemble from the kinetically accessible states of the 
unbound folded protein. Consequently, the kinetics of large scale conformational transi- 
tions between two meta-stable states are determined largely by the inherent conformational 
dynamics within the folded state free energy basin. Such conformational dynamics imply 
an inherent flexibility or "intrinsic plasticity" of the folded state. In this paper, we focus 
on how this inherent flexibility of a protein molecule influences the mechanism controlling 
the kinetics of allosteric transitions. There are a couple of possible scenarios for the tran- 
sition mechanism in terms of changes in conformational flexibility. One possibility is that 
the flexibility adjusts smoothly to the conformation deformation between two specific meta- 
stable states in the free energy surface. In this case, the inherent flexibility of the protein 
remains relatively constant if the meta-stable states have similar flexibilities or else changes 
smoothly between the flexibilities of the meta-stable states if the flexibilities are different. 
Another possible mechanism, called "cracking",- has recently been proposed as an alter- 
native mechanism that may be important in some conformational transitions. In terms of 
conformational flexibility, cracking involves non-monotonic changes in flexibility along the 
transitions route. In particular, the flexibility of specific regions of the protein may tran- 
siently increase through local unfolding. As explained in RefJ2|, local unfolding can relieve 
specific areas of high stress during conformational deformations that would result in high 
free energy barriers if the protein remained uniformly folded throughout the transition. 

In the model developed to explore this idea, cracking is introduced directly into the 
formalism as means to incorporate nonlinear elasticity in an otherwise harmonic description 
of conformational fluctuations (i.e., normal modes) .-"^ This work clearly shows that local 
unfolding can dramatically lower predicted free energy barriers between local free energy 
minima and hence facilitates faster kinetics. In contrast, the variational model developed in 
the present work does not assume cracking from the outset. The evidence of cracking here 
arises entirely from the analysis of the local changes in flexibility along transition routes 
predicted by an inherently non-linear model of conformational transitions. Our results are 



thus an independent verification of tlie ideas and formalism developed in the model presented 
in Ref yy. 

Cracking does not occur in all conformational transitions, but even when it does, the 
average conformations of the local minima in the free energy may show little signs that 
local unfolding is involved.- Instead, cracking is a subtle consequence of the nature of the 
conformational deformation required to connect the two states. Since cracking is akin to 
folding of the whole protein albeit on a constrained and local scale, it is reasonable to ask 
if some insight appropriated from the success of the energy landscape theory of protein 
folding^ can help to anticipate when cracking is likely to be important in conformational 
transitions between two distinct folded conformations. One motivation for the work in this 
paper is to investigate if the cracking mechanism of conformational transitions, like the 
mechanism for folding of two state proteins,- is determined by structural topology of the 
two meta-stable states. To address this question we study the conformational transitions of 
the open and closed conformations of the two homologous domains of calmodulin (CaM) as 
well as a fragment of CaM involving parts of both domains. Aside from general theoretical 
implications, the results presented in this paper are interesting from the point of view of 
understanding the well-known and intensely studied thermodynamic differences between the 
two domains of CaM.-i^ 

Our approach is based on a coarse-grained variational model previously developed to 
characterize protein folding.^ This model is in harmony with several recent coarse-grained 
simulations based on a folded-state biased interaction potential that interpolates between 
the contact maps of the meta-stable states.-'^'^'^^'i^i^ Interestingly some of these recently 
proposed coarse-grained model for protein conformational transition can also capture the 
cracking or local partial unfolding during transition^^ii^ii^. In this paper we compare the 
detailed predicted transition routes for the open/closed conformational change of the C- 
terminal domain of Calmodulin (cCaM) with that of the N-terminal domain of CaM (nCaM). 
The open and closed state of each domain is shown in Fig. [T^a-b). We find that although the 
nCaM and cCaM have similar structures, the calculated transition routes predict that the 
cracking is essential only in cCaM. It is well established that the two domains differ in flexi- 
bility as well as binding affinity for Ca^^.^'^ These properties are not determined by topology 
(since they are the same) but ultimately by more subtle structural differences encoded in 
the sequence of the two domains. Our results suggest that nCaM is flexible enough that 



cracking is not necessary to relieve stress during the conformational change, while the cCaM 
is relatively rigid so that cracking is involved in essentially the same conformational change. 
In addition, we investigate the conflicting Ca^+ binding mechanism in the "odd-even" and 
"even-odd" paired EF-hands between the two domains of CaM [shown in Fig. [11(c)]. Our 
study indicates that the highly flexible central linker region locally unfolds before folding 
into a stable a-helix. Our predicted transition mechanism rationalizes the sequential Ca^^ 
binding in CaM2/3.— in terms of the differences in conformational flexibility of Ca^+-binding 
loops in CaM2/3. Finally, we address how our comparative study of conformational transi- 
tions in "odd-even" and "even-odd" paired EF-hands reveals the interplay between intrinsic 
plasticity, target binding affinity, and function of CaM. 

RESULTS 

Conformational Flexibility of the CaM Domains. We first compare the inherent 
flexibility of the two domains of CaM by calculating the mean-square fluctuations Bi for 
each residue in the nCaM and cCaM for the open/closed conformational transitions of each 
domain. Fig. [2](a) shows Bi along the open/closed conformational transition for each residue 
of nCaM from our very recent studies.— The magnitude of the fluctuations for each residue 
of cCaM for the open/closed transition route is shown in Fig. |3](a). Comparison of Fig. [2|^a) 
and [21(a) show that although the binding loops and helix-linker in both domains are very 
flexible, the apo (closed)-nCaM is inherently more flexible than apo (closed)-cCaM. For 
nCaM, the flexibility of the the binding loops I and II and the B/C helix-linker between EF- 
hands 1 and 2 decreases upon domain opening. Similarly, the flexibility of the Ca^+ binding 
loops III and IV decreases considerably during the opening of cCaM. The fluctuations of the 
F/G helix-linker between EF-hands 3 and 4 has different behavior along the transition route. 
The flexibility of this region of the protein increases to a relatively high flexibility in the 
intermediate states between the open and closed conformations before reducing the flexibility 
of the folded open conformation. (Fig. [3]^a)). We also note that the F/G helix-linker in the 
open structure is more flexible than in the closed structure. The calculated conformational 
flexibility suggests that unlike the B/C helix-linker in the apo (closed)-nCaM the F/G helix- 
linker in the apo (closed)-cCaM is relatively less flexible. The increase and then decrease 
in flexibility of the F/G helix-linker during the domain opening of cCaM indicates there is 



cracking or local partial unfolding in this less flexible F/G helix-linker which facilitates the 
conformational change by increasing its flexibility near the transition state ensemble. [See 
the change in fluctuations of the F/G helix-linker during domain opening in Fig. [21(a)] The 
model also predicts that binding loop II in CaM has the highest flexibility in comparison 
with other binding loops. [Compare the flexibility of the binding loops in Fig. [2]^a) and 
[3]^a).] This result agrees with a recent molecular dynamics simulation study of CaM D129N 
mutant, where the binding loop II consistently exhibited higher mobility than the other 
three Ca^"*" binding loops.— 

Cracking in the Conformational Transition of cCaM. To compare the confor- 
mational flexibility of the two CaM domains in more detail we have plotted the change 
in fluctuations Bi for some speciflc residues in Fig. HJ Each binding loop of CaM has 12 
residues with a conserved Glutamate (Glu) residue in the 4^^ position of the loops. In Fig. H] 
we show Glu23 and Glu59 from binding loops I and II of nCaM, respectively. Also, Glu96 
and Glul32 from binding loops III and IV of cCaM, respectively. In Fig. |H we also included 
the Glycine (Gly) and Aspartic (Asp) residues from the helix-linker at position 45 and 118 
of the nCaM and cCaM, respectively. [Note that all these six residues show sharp peaks in 
Fig.[2](a) and [3t^a).] Residue AspllS of the F/G helix-linker from cCaM shows very different 
behavior relative to the rest of the residues in Fig. |H Residue AspllS shows highest Bi near 
the transition state at ao = 0.4 whereas Bi of all the other residues shown decrease mono- 
tonically during closed to open transition with the closed conformation having the largest 
fluctuations. This result supports that this increase in flexibility of the helix-linker region 
in cCaM near the transition state is due to some local transient unfolding or cracking. For 
further analysis of cracking in the cCaM we have also compared the contact pair potential 
energy Uij for residues Glu45 and AspllS from the helix-linker of nCaM and cCaM respec- 
tively along the transition route. Fig. [5] shows that for residue AspllS, Uij increases in the 
transition state and therefore cracking occurs. On the other hand, Uij for residue Glu45 
from the helix-linker of nCaM deceases in the transition state. This striking difference in 
Uij implies that the increase in contact energy of the F/G helix-linker during opening of 
cCaM leads to some transient cracking in this region near the transition state which further 
enhances the inherent flexibility of this linker during the open/closed conformational change. 

Open/Closed Transition Mechanism of the CaM Domains. To further elucidate 
the predicted conformational transition mechanisms of domain opening in CaM, we consider 



a structural order parameter that measures the similarity to the open (holo) or closed (apo) 
state conformations, Api', described in the Methods section. This order parameter is defined 
such that Api" = 1 corresponds to the closed conformation and A'pl = — 1 corresponds to 
the open conformation of nCaM or cCaM. The transition routes illustrated by Api for each 
residue of nCaM and cCam are shown in Fig. El^b) and Fig. [3](b), respectively. Fig. EJ^b) 
illustrates the conformational transition of cCaM in terms of Ap^ for each residue. As 
shown in Fig. [3](b), the model predicts that the structural change in binding loop II occurs 
earlier than binding loop III during the domain opening of cCaM. [See the sharp change 
in color near residue number 130 in Fig. [3](b)]. Also, helix G has an earlier conformational 
transition than helix F (Fig. El^b)). This sequence is quite different from nCaM as shown 
in Fig. [2](b). For nCaM the structural change in helix B is predicted to be earlier than the 
structural change in helix C for the closed to open conformational transition. Finally, when 
we compared the structural change in the helix-linker of the two domains, our results show 
F/G helix-linker in the cCaM has abrupt transition near the open state. [See the transition 
near residue number 115 in Fig. [3](b).] This may have some implication for cracking in the 
F/G helix-linker region of cCaM. In contrast, the transition in the B/C helix-linker in nCaM 
is more gradual (see Fig. [2](b)). 

Conformational Flexibility and Cracking of the CaM2/3 Fragment. The en- 
gineered "even-odd" EF-hands paired CaM2/3 fragment (composed of EF-hands 2 and 3) 
has been shown recently to have distinct transition characteristics from the "odd-even" 
EF-hands paired nCaM and cCaM.— NMR spectroscopy has shown that Ca^+-free (apo) 
CaM2/3 does not have a stable folded structure but rather shows characteristics of a molten 
globule state. Ca^"*" binding induces the folding of this "even-odd" paired EF-hand motifs 
CaM2/3 and the Ca^"'"-bound (holo) CaM2/3 adopts a similar structure as holo-nCaM or 
cCaM (see PDB ID 2hf5)M In contrast to the Ca2+-CaM (holo-CaM) structure, the Ca2+- 
CaM2/3 structure does not have a stable helix in the central linker region. In our model, 
we study the two meta-stable conformations, apo-CaM2/3 and holo-CaM2/3 were taken 
directly from the folded apo-CaM (PDB ID Icfd) and holo-CaM (PDB ID Icll) respectively. 
Here, we focus mainly on the conformational flexibility and transition of the central linker to 
a stable a-helix. This central linker of CaM was also previously studied by MD simulation.— 
Even though our model does not accommodate all the relevant conditions of the NMR ex- 



periment reported in Ref. Il5|, the model does capture certain aspects of the the apo/holo 



conformational transition of this CaM2/3 fragment. In particular, the model predicts a apo- 
CaM2/3 to the holo-CaM2/3 conformational transition mechanism that agrees well with the 
sequential Ca^^ binding mechanism of CaM2/3 suggested by the NMR measurements. 

The magnitude of the fluctuations of each residue in CaM2/3 [shown in Fig. [Hl^a)] suggests 
that the change in conformational flexibility of binding loops II and III of CaM2/3 are 
very different. In particular, helix C and binding loop II are highly flexible in the apo 
conformation of CaM2/3. This very high intrinsic flexibility of apo-CaM2/3 may account 
for the molten globule state characteristics described in Ref. Il5|. Fig. [6](a) also shows that 
binding loop II is more flexible than loop III. This result is in harmony with the NMR 
measurements of Ca^^ binding in CaM2/3^ that shows Ca^^ binds to loop III with higher 
affinity (lower Ca^+ concentration) than loop II. Our results also indicate that the flexibility 
of binding loop III decreases monotonically along the transition route from the apo-CaM2/3 
to the holo-CaM2/3 structure. In contrast, the flexibility of helix F increases along the 
transition route from apo to holo transition of CaM2/3 fragment. In particular, N-terminal 
part (close to binding loop III) of helix F partially unfolds when adopting the holo-CaM2/3 
conformation. We also notice in Fig. [6](a) that the domain linker between EF-hands 2 and 
3 (residue number 74-81) during the apo to holo transition has a relatively large change 
in conformational flexibility. The fluctuation amplitude Bi of this linker increases in the 
intermediate states and then decrease abruptly as the helix forms near the holo state. The 
local unfolding of the domain linker enhances the flexibility of this linker region dramatically, 
first relaxing in structure and then stabilized in the holo-CaM2/3 structure as a rigid a-helix 
conformation. This local unfolding signaled by increase and decrease in flexibility is similar 
to the cracking exhibited in cCaM. 

Local intermediate unfolding (i.e., cracking ) of the domain linker during the apo to 
holo transition of CaM2/3 can also be seen through change in the pair potential energy 
Uij for the Aspartic (Asp) residue from the linker at position 80 along the transition route 
as shown in Fig. [71 During the apo to holo transition of CaM2/3 the energy of residue 
Asp80 increases from the apo-CaM2/3 structure {a = 1) and decreases abruptly at the 
transition state (ao = 0.4) to the minimum at holo-CaM2/3 structure (ao = 0) as the helix 
forms . This increase and then sharp decrease in the contact energy of residue Asp80 signals 
cracking (local unfolding and refolding) of the linker in CaM2/3 to a stable a-helix. We note 
also that the energy relative to to the holo-CaM2/3 energy of Asp80 is higher than that of 



AspllS from the F/G helix-linker of cCaM. This emphasizes the importance of cracking in 
the apo-CaM2/3 to holo-CaM2/3 transition. 

Conformational Transition Mechanism of CaM2/3 Fragment. For CaM2/3, the 
order parameter Ap7 = 1 corresponds to the apo state and Ap7 = — 1 corresponds to the holo 
state of CaM2/3. The abrupt and early change in A'pl = 1 for residues in the domain linker 
(sequence number 74-81) illustrated in Fig. [6](b) clearly shows the transition of the flexible 
inter domain linker to a rigid a-helix as discussed above in terms of the fluctuations {Bi}. 
Fig. Mjo) also indicates that the structural change in binding loop III is initiated earlier than 
the structural change in binding loop II. Also, the conformational change in binding loop 
II is much more gradual than that of binding loop III. Similar to the identification of the 
Ca^^ binding affinity with the changes in inherent flexibility, the early structural change in 
binding loop III may imply the stepwise binding of Ca^"^ in CaM2/3 with loop III having a 
higher binding affinity than loop II. 

DISCUSSION 

CaM, a small (148- residue) Ca^+-binding protein with very high plasticity, may be an 
ideal system to demonstrate flexibility influenced conformational transitions. The protein 
consists of structurally similar N- and C-terminal globular domains connected by a flexible 
tether also known as central or interdomain linker.— The Ca^"*" induced structural rear- 
rangement in CaM result in the solvent exposure of large hydrophobic surface responsible 
for molecular recognition of various cellular targets.— 

While similar in structure and fold, the two CaM domains are quite different in terms of 
their flexibility, melting temperatures, and Ca^^-binding affinities. -^'^ In the two homologous 
(46% sequence identity) domains of CaM, Ca^"*" binding occurs sequentially. First, in the 
binding sites of cCaM and then in the binding sites of nCaM.- In spite of large cooperativ- 
ity in the Ca^"*" binding process within each domain the two domains reflect different Ca^+ 
and target affinities. The N-terminal pair of EF-hands binds to Ca^+ ions with much lower 
affinity than the C-terminal EF-hands pair.- Nuclear Magnetic Resonance (NMR)^'2i and 
molecular-dynamics (MD) simulation^^ studies have shown that the Ca^^-bound nCaM is 
considerably less open than the cCaM. This was not observed in the X-ray crystal structure 
of the protein. Experimentally it has been shown that the more-conserved cCaM has a 
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greater affinity for Ca^"*" and some CaM targets, whereas the nCaM is less specific in its 
choice of target motif.— 12^ Heat denaturation studies have shown that cCaM of Ca^'^-free 
(apo) CaM starts to denature shghtly above the physiological temperature.^ The denatura- 
tion of the apo-cCaM was observed at lower concentration of denaturant than denaturation 
of the nCaM, while the order was reversed for Ca^+-CaM.^^ From temperature-jump fluo- 
rescence spectroscopy by Rabl et a\M the instability of the cCaM was also observed from 
the study of unfolding of apo-CaM. They suggested that the cCaM was partially unfolded 
at native conditions. Recent NMR experiments done by Lundstrom and Akke^^ monitoring 
relaxation rates involving ^^C" spins in adjacent residues of E140Q mutant of cCaM revealed 
transient partial unfolding of helix F. This was interpreted as a global exchange process in- 
volves a partially unfolded minor state that was not detected previously2^i2E. A very recent 
conformational dynamics simulation of the cCaM by Chen et al.— has shown presence of an 
unfolded apo-state. These observations further suggests that the conformational exchange 
may be more complex than a simple two-state process. 

The variational model presented in this paper agrees with the emphasis on differences in 
inherent flexibility of the two domains to understand their distinct physical characteristics 
and the complexity in the conformational transition mechanism.-^ Our results also agree 
with a MD simulation study'^^ which has shown that the nCaM is inherently more flexible 
with lower binding affinity than the cCaM. The more open conformation and lower intrinsic 
flexibility of the cCaM is also probably the key to understanding initial binding between 
this domain and CaM's target enzymes.—'^ 

The variational model presented in this paper predicts that the different inherent flexi- 
bilities of the two domains of CaM also lead to distinct transition mechanisms, even though 
the folded state topology of the two domains is the same. The mechanism controlling the 
open/closed transition does not involve cracking for the relatively flexible nCaM whereas the 
mechanism controlling the conformation transition of the more rigid cCaM exhibits tran- 
sient partial local unfolding in its helix-linker between the binding loops. This partial local 
unfolding or cracking in the cCaM shows the complexity of the open/closed conformational 
transition mechanism of cCaM. 

A recent NMR experiment studied the EF-hand association affects on the structure, 
Ca^+ affinity, and cooperativity of CaM.^^ The EF-hands in CaM domains are "odd-even" 
paired with EF-hands 1 and 2 in nCaM and EF-hands 3 and 4 in cCaM. This arrangement 



of EF-hands in CaM is thought to be a consequence of its evolution from a biologically 
related ancestor EF-hand by gene duplication.— The "even-odd" pairing of EF-hands 2 and 



3 (CaM2/3) (Fig. [I](c)) has been characterized by NMR in Ref. [ISJ. In this fragment, EF- 
hands 2 and 3 are connected by the central linker between nCaM and cCaM. Although, from 
the crystal structure of holo Ca^+-CaM this interdomain linker region is observed to be a 
long rigid a-helix,— several NMR relaxation experiments have demonstrated that this central 
linker is flexible in solution near its midpoint and the two domains do not interact.— ^^^ In 
contrast to the high affinity and positive cooperativity for Ca^"*" binding in the two "odd- 
even" paired EF-hand domains of CaM, the CaM2/3 binds Ca^"*" in a sequence. First in 
the high-affinity EF-hand 3 and then in the EF-hand 2 with much lower affinity. ^^ Although 
not a direct focus of our paper, we also note that a peptide binding to CaM2/3 has also 
been characterized recently by Mcintosh and co-workers.^^^ It was found that CaM2/3 adopts 
Ca^+-bound structure with peptide binding very much similar to those of nCaM or cCaM. 
These observations reflect the very high plasticity of the EF-hand association and mediate 
Ca^+-dependent recognition of target proteins. 

Our results of apo/holo conformational transition of CaM2/3 reveals that the central 
linker in CaM2/3 is highly flexible and dominates the transition mechanism. The folding of 
this linker to a rigid a-helix is preceded by a further increase in its flexibility and local un- 
folding. This cracking leads to a very sharp transition of the helix-linker from apo-CaM2/3 
to holo-CaM2/3 conformation. While our model did not predict partial local unfolding or 
cracking of helix-F from the open/closed conformational transition of cCaM, the apo/holo 
conformational transition of CaM2/3 reveals cracking in some region of helix-F as we have 
discussed already in comparison with the NMR study^. The predicted change in conforma- 
tional flexibility of CaM2/3 reveals the stepwise binding of Ca^"*" ions with binding loop III 
having higher affinity than binding loop II. Furthermore, our results suggests that due to 
very high inherent flexibility the central linker in CaM2/3 behaves more like a helix-linker in 
CaM domains, as concluded from NMR measurements.— On the other hand, in Ca^^'^-bound 
(holo) CaM structure^^ this same central linker adopts a more stable a-helix conformation 
which requires cracking in the conformational transition from its apo-CaM structure. 
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METHODS 

Preparation of Structures. The structures of the two conformations of the transition 
are rotated to have the same center of mass and minimum root-mean square deviation of the 
Ca positions. In nCaM, EF-hands 1 and 2 are consist of hehces A/B and C/D (Fig. [11(a)), 
respectively, whereas hehces E/F and G/H in cCaM (Fig. [11(b)) form EF-hands 3 and 4, 
respectively. The fragment CaM2/3 contains EF-hands 2 and 3 with helices C/D and E/F 
(Fig.[Tl(c)), respectively. The conformational transition in cCaM is modeled from the residues 
76-147 of apo-CaM (PDB ID Icfd)^ and holo-CaM (PDB ID IcU) ^ structures, while the 
conformational transition in nCaM is modeled from the residues 4-75 of the same PDB 
structures. The apo to holo conformational transition of the "even-odd" paired EF-hand 
motifs CaM2/3 are modeled from EF-hands 2 and 3, residues 46-113 of apo-CaM (PDB ID 
Icfd) and holo-CaM (PDB ID Icll) structures. 

The Variational Model of Conformational Transitions. A conformation of a 



protein in our model (described in more detail in Refll6l) is described by the A^ position 
vectors of the a-carbons of the polypepetide backbone, {tj}. Partially ordered ensembles of 
polymer configurations are described by a coarse-grained reference Hamiltonian 

Ho/kBT = H^^^Jk^T + A ^ Q[r, - rf (a,)]^ (1) 

i 

where T is the temperature and k-Q is Boltzmann's constant. Here, Tichain is a harmonic 
potential that enforces chain connectivity of a freely rotating polymer with mean bond 
length a = 3.8A,— and mean valance angle cos 6* = 0.8.'^^ The second term includes N 
variational parameters, {Cj}, that control the magnitude of the fiuctuations about a-carbon 
position vectors 

rf(a.) = a.rf- + (l-«,)rf-°. (2) 

Here, we have introduced another set of A^ variational parameters, {ai} ( ranging between 
and 1), that specify the backbone positions as an interpolation between the apo-CaM and 
holo-CaM conformations, {r^ ''^°} and {r^ hoio j^ respectively. The probability for a particular 
configurational ensemble specified by the variational parameters {Cj, a^} at temperature T 
is given by the variational free energy F{{C}, {a}) = E{{C}, {a}) — TS{{C}, {a}). Here, 
S{{C}, {a}) is the entropy loss due to the localization of the residues around the mean 
positions {ri{ai)}. The energy is determined by the two-body interactions between distant 
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residues E{{C},{a}) = ^uj]^ijUij, where Uij is the average of the pair potential u{rij) 
over Ho,— and e^j is the strength of a fully formed contact between residues i and j given 
by Miyazawa-Jernigan interaction parameters.— The sum is restricted to a set of contacts 
[i,j] determined by pairs of residues in the proximity in each of the meta-stable confor- 
mations.-- Each meta-stable conformation has a distinct but overlapping set of contacts 
[hJ]a.po, [^, j]hoio and [i,j] is the union of these sets of contacts. In this model, the inter- 
action energy for contacts that occur exclusively in only one meta-stable structure is given 
by exp(— Mij/Z^B^) = 1 + exp(— (-Uij^/A^B^)- This two-state model provides an interpolation 
at the individual contact level of the interaction energy determined by the two meta-stable 
conformations and is similar to some other native state biased potentials describing confor- 
mational transitions.— 'ii'i^ 

Analysis of the free energy surface parameterized by {C, a} follows the program developed 
to describe folding:— the mechanism controlling the kinetics of the transitions is determined 
by the ensemble of structures characterized by the monomer density at the saddlepoints of 
the free energy. At this point, we simplify our model and restrict the interpolation parameter 
ftj to be the same for all residues, «« = ao-~ With this simplification, the numerical problem 
of finding saddle-points with respect to {C, a} simplifies to minimizing the free energy 
F{{C}, ao) with respect to {C} for a fixed oq. 

Conformational Flexibility. Main-chain conformational flexibility is characterized by 
the mean-square fluctuations {Bi{ao)} of each a-carbon of the polypepetide chain from its 
average positions, 6ri, along the transition route as a goes from to 1.— These natural order 
parameters for the reference Hamiltonian Tio, Bi = (5r^)o, contains information about the 
degree of structural order of each residue.-"^ 

Conformational Transition Mechanism. The main-chain dynamics responsible for 
the detailed mechanism of the conformational transition in our analysis is based on the 
order parameters introduced to describe the folding mechanism p^^° °° = (exp(— a'^(rj — 
rj ='p°/''^°'°)2^^ with a^ = 0.5 defining the width of a Gaussian window about the meta-stable 
structure {r^ apo/hoio|^ j^ particular, it is convenient to characterize the relative similarity to 
the apo structure along the transition route through the normalized measure 

-apo/ N _ Pr°(«o)-Pr°(0) /o^ 

^' 1«0J ^apo^^^_^apo^Q^ , {<i) 

where p^'^°(q;o) is the monomer density of the i^^ residue with respect to the apo conformation 
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described by {r^ ^p°yi6_ Similarly, we represent the relative structural similarity to the holo 
conformation as 

^holo. ^ ^ Pr°K)-Pr°(l) U) 

P^ 1«0J ^holo(o)_pholo(i) ' W 

where p^°'°(ao) is the monomer density of the i^^ residue with respect to the holo conforma- 
tion described by {r^ hoio j_ j^^ ^^^e holo state, pi''^P°(0) = and Pi'^°'°(0) = 1, while in the apo 
state /^^^^"(l) = 1 and 'pi^°^°{l) = 0. To represent the structural changes more clearly, it is 
convenient to consider the difference, 

AA(ao)=r''°(«o)-pl'°'°M (5) 

for each residue. This difference shifts the relative degree of localization to be between 
Ap^(l) = 1 and Ap7(0) = — 1 corresponding to the apo and holo conformations, respectively. 
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Figure Legends 

Figure {Jl 

Three dimensional structures of calmodulin (CaM) domains and EF-hands 2 and 3 frag- 
ment. The apo-CaM and holo-CaM structures shown here are correspond to human CaM 
with PDB code Icfd and ic// respectively, (a) The closed, apo and open, holo conformations 
of N-terminal domains of CaM (nCaM) consist of helices A/B and C/D with binding loops I 
and II respectively, (b) The closed, apo and open, holo conformations of C-terminal domains 
of CaM (cCaM) consist of helices E/F and G/H with binding loops III and IV respectively, 
(c) The unfolded (apo) and folded (holo) EF-hands of CaM (CaM2/3) consist of helices C/D 
and E/F with binding loops II and III respectively. These three-dimensional illustrations 
were made using Visual Molecular Dynamics (VMD).— 

Figure [3 

(a) Fluctuations Bi vs residue index of N-domain of CaM (nCaM) for selected values 
of the interpolation parameter ao in the conformational transition route between open and 
closed structures. The secondary structure of nCaM is indicated above the plot. Helices are 
represented by the rectangular boxes, binding loops and helix-linker are by lines and small 
/^-sheets are by arrows, (b) Difference between the normalized native density Ap^ (a measure 
of structural similarity) of each residue for different Oq- The change in color from red to blue 
is showing the closed to open conformational transition of nCaM. This is normalized to be 
— 1 at the open state minimum (oq = 0; blue) and 1 at the closed state minimum (oq = 1; 
red). 

Figure O 

(a) Fluctuations Bi vs residue index of C-domain of CaM (cCaM) for selected values of the 
interpolation parameter ao in the conformational transition route between open and closed 
structures. The secondary structure of cCaM is indicated above the plot, (b) Difference 
between the normalized native density Api (a measure of structural similarity) of each 
residue for different a^. The change in color from red to blue is showing the closed to open 

16 



conformational transition of cCaM. This is normalized to be —1 at the open state minimum 
(ao = 0; blue) and 1 at the closed state minimum (ao = 1; red). 

Figure [^ 

Change in fluctuations Bi of the residues from binding loops and helix-linker of CaM 
domains along the open/closed conformational transition route for different Oq. Residues 
Gly23, Gly59 and Glu45 are from binding loop I, loop II and the helix-linker of the nCaM, 
respectively. While, residues Gly96, Glyl32 and AspllS are from binding loop III, loop IV 
and the helix-linker of the cCaM, respectively. 

Figure [3 

Average pair potentials Uij of the residues from the helix-linker of CaM domains along 
the open/closed conformational transition route for different a^. Residue Glu45 and AspllS 
are from the nCaM and cCaM, respectively. 

Figure O 

(a) Fluctuations Bi vs residue index of fragment EF-hands 2 and 3 of CaM (CaM2/3) 
for selected values of the interpolation parameter a^ in the conformational transition route 
between unfolded and folded states. The secondary structure of CaM2/3 is indicated above 
the plot, (b) Difference between the normalized native density Ap7 (a measure of structural 
similarity) of each residue for different a^. The change in color from red to blue is showing 
the unfolded to folded conformational transition of CaM2/3. This is normalized to be —1 
at the folded state minimum (ao = 0; blue) and 1 at the unfolded state minimum (ao = 1; 
red). 

Figure [^ 

Average pair potentials Uij of the residues AspSO from the central linker of CaM along 
the open/closed conformational transition route for different ao- 
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